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ABSTRACT 



We determine the yields of the elements from Na to Ni for Type II supernovae 
(SNe II) and the yield patterns of the same elements for Type la supernovae (SNe 
la) and very massive (> 100 M Q ) stars (VMS) using a phenomenological model of 
stellar nucleosynthesis and the data on a number of stars with —4 < [Fe/H] < —3, 
a single star with [Fe/H] = —2.04, and the sun. We consider that there are two 
distinct kinds of SNe II: the high-frequency SNe ll{H) and the low-frequency 
SNe H(L). We also consider that VMS were the dominant first-generation stars 
formed from big bang debris. The yield patterns of Na to Ni for SNe 11(H), ll(L), 
and la and VMS appear to be well defined. It is found that SNe 11(H) produce 
almost none of these elements, that SNe II (L) can account for the entire solar 
inventory of Na, Mg, Si, Ca, Ti, and V, and that compared with SNe ll(L), VMS 
underproduce Na, Al, V, Cr, and Mn, overproduce Co, but otherwise have an 
almost identical yield pattern. A comparison is made between the yield patterns 
determined here from the observational data and those from ab initio models of 
nucleosynthesis in SNe II and VMS. 

The evolution of the other elements relative to Fe is shown to involve three 
distinct stages. The earliest stage is in the domain of [Fe/H] < —3 and is gov- 
erned by VMS activities with some small contributions from SNe II, all of which 
are dispersed in a dilution mass of ~ 10 6 -10 7 M Q . The beginning of the 

second stage is marked by the cessation of VMS activities and the onset of major 
formation of normal stars (with masses of ~ 1-60 M ) at [Fe/H] —3. The dilu- 
tion mass for SN II contributions also drops sharply to Mj^ 11 « 3 x 10 4 M & . The 
subsequent quasi-continuous chemical evolution until [Fe/H] ~ —1 is governed 
by SNe 11(H), which produce mainly the heavy r-process elements above Ba, 
and SNe H(L), which produce essentially all the other elements. The third stage 
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starts with the onset of SN la contributions to mainly the Fe group elements at 
[Fe/H] ~ —1. The domain of [Fe/H] > —1 is then governed by contributions 
from SNe II(-P), H(L), and la and other low mass stars. It is shown that the 
abundances of non-neutron-capture elements in stars with [Fe/H] < and those 
of r-process elements in stars with [Fe/H] < — 1 can be well represented by the 
sum of the distinct components in the phenomenological model. The proposed 
evolutionary sequence is directly related to the problems of early aggregation 
and dispersion of baryonic matter and to the onset of formation and chemical 
evolution of galaxies. It is argued that the prompt inventory governed by VMS 
contributions should represent the typical composition of dispersed baryonic mat- 
ter in the universe. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — Galaxy: evo- 
lution — Stars: Population II 

1. Introduction 

In this paper we discuss the nucleosynthetic yields of supernovae (SNe) and very massive 
stars (VMS, with masses of > 100 M Q ) based on the observed abundances in stars with 
—4 < [Fe/H] < — 1. It was considered in our previous works [Wasserburg & Qian 2000a 
(WQ),b; Qian & Wasserburg 2001a,c (QW)] that VMS were formed from big bang debris and 
provided an initial or prompt (P) inventory of Fe and associated elements up to [Fe/H] —3. 
It was assumed that normal stars (with masses of ~ 1-60 M Q ) could only be formed at 
[Fe/H] > —3. Subsequent evolution of normal stars led to Type II SNe (SNe II) and at 
much later times (corresponding to [Fe/H] > —1) also to Type la SNe (SNe la), all of 
which provided further chemical enrichment of the interstellar medium (ISM) beyond the P- 
inventory. The special status of [Fe/H] m —3 in chemical evolution of the ISM was identified 
based on a sharp increase in the abundances of heavy r-process elements (Ba and above) at 
this metallicity (WQ). This is shown in Figure 1 for Ba. The sharp increase was attributed to 
the rapid occurrence of high-frequency SNe II(P) that produce mainly the heavy r-process 
elements but no Fe. The Fe enrichment of the ISM at —3 < [Fe/H] < — 1 was attributed 
to the low-frequency SNe II (L) and that at [Fe/H] > -1 to both SNe II (L) and la. In 
addition to contributing 1/3 of the solar Fe inventory, SNe H(L) are mainly responsible 
for the light r-process elements (Ba and below). QW have derived the P-inventory and SN 
II(P) and ll(L) yields of r-process elements from the observed abundances in two stars with 
[Fe/H] ~ —3 and the solar r-process abundances. They have shown that the abundances 
calculated from the three- component model including the P-inventory and the contributions 
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from SNe 11(H) and ll(L) were in good agreement with several independent data sets for 
a large number of r- process elements over the wide range of —3 < [Fe/H] < —1. The 
same agreement was also achieved for Sr, Y, Zr, and Ba when the standard solar r-process 
abundances of these elements (Arlandini et al. 1999) were increased by factors of ~ 2-6. In 
general, the model of QW gave a satisfactory description of the chemical evolution of all the 
r-process elements over —3 < [Fe/H] < —1, which corresponds to a period of ~ 3 x 10 9 yr 
in the early Galactic history. 

Here we extend the model of QW to the elements from Na to Ni and derive the relevant 
yields of SNe II(-H') and ll(L) and the yield pattern of SNe la (§2). We also infer the yield 
pattern of VMS from the observed abundances in stars with —4 < [Fe/H] < —3 (§3). It will 
be shown that SNe 11(H) produce almost none of the elements from Na to Ni, that SNe H(L) 
can account for the entire solar inventory of Na, Mg, Si, Ca, Ti, and V, and that compared 
with SNe II (L), VMS underproduce Na, Al, V, Cr, and Mn, overproduce Co, but otherwise 
have an almost identical yield pattern. We re- interpret the metallicity of [Fe/H] —3 as 
corresponding to the cessation of VMS activities and the onset of major formation of normal 
stars. It will be shown that VMS dominated the production of Na to Ni at [Fe/H] < —3 
and that explosions of these objects could induce mixing within ~ 10 6 -10 7 M & of ISM, to 
be compared with «3x 10 4 M q for SNe II. A model for general chemical evolution of the 
ISM starting from a zero metallicity state will be presented (§4). The wide range in yields 
calculated from nucleosynthetic models for SNe II and VMS will be discussed in light of 
the rather constant and metallicity-independent yield patterns as determined here from the 
observed abundances in metal-poor stars (§5). The cosmological implications of the results 
presented here will also be discussed (§5). 



2. Three-Component Model and the P-Inventory and SN Yields 

We consider a homogeneous mass of gas and assume that stars formed from this gas have 
the composition of the gas. The chemical evolution of an element E in this gas is specified by 
the initial number of E atoms per H atom in the gas and the number of E atoms per H atom 
added from various SNe II to the gas. As shown by QW, if each distinct kind of SN II has a 
fixed relative yield pattern, then the inventory of E atoms may be simply described by the 
effective numbers of SNe II that contributed to the gas. The effective numbers are calculated 
by assuming constant yields and a standard dilution mass for SNe II. The results of astration, 
fragmentation, and merging of different gas masses, as well as the effects of variable SN II 
yields (but with fixed yield patterns), are subsumed in the representation of the effective 
numbers of contributing SNe II. By using this approach, it was shown that the abundances 
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of r-process and related elements (Sr and above) in stars with —3 < [Fe/H] < —1 can be 
described by a three-component model including the P-inventory and the contributions from 
SNe 11(H) and 11(f) (QW). For a star formed from an ISM with a number np of contributing 
SNe II (if) and a number Ul of contributing SNe 11(f), the abundance of an element E in 
the star is given by 

10 log6(E) = 10 lo gep (E) +UhX 1Q log6 ff (E) + Ul x 10 lo geL (E) ? ^ 

where logep(E), logep(E), and logeL(E) represent the P-inventory and the SN 11(H) and 
11(f) yields of E, respectively. Here and below, the standard spectroscopic notation log e(E) = 
loge(E/H) + 12 is used. The parameter logep(E) or loge£,(E) corresponds to the abundance 
of E in the ISM resulting from a single SN 11(H) or 11(f) for a standard dilution mass 
(Mf^ 11 3 x 1O 4 M ) of baryonic matter that is free of "metals" and dominantly composed 
of H. This dilution mass is fixed by the frequency of SNe II(ff ) as required for replenishment 
of fresh radioactive 182 Hf in the ISM (e.g., Wasserburg, Busso, & Galllino 1996; QW) and 
is consistent with the typical total amount of ISM swept up by an SN II remnant (e.g., 
Thornton et al. 1998). 

As the heavy r-process elements above Ba, e.g., Eu, are exclusively produced by SNe 
II(ff), the number n H for the star can be obtained from its observed loge(Eu): 

Uh _ ^gloge(Eu)-loge H (Eu) _ -j^glog e(Eu)+2.48 ^ ^) 

where logep(Eu) = —2.48 (QW) is used. Likewise, as additions of Fe at —3 < [Fe/H] < —1 
beyond the P-inventory are made only by SNe 11(f), the number til for the star can be 
obtained from its observed loge(Fe): 

n L = [io lo ^( Fc ) - io lo s^( Fc )] /io lo s^( Fe ), (3) 

where loge P (Fe) = 4.51 and loge L (Fe) = 5.03 (QW). By using [Fe/H] = loge(Fe)-loge (Fe) 
with loge (Fe) = 7.51 (Anders & Grevesse 1989), equation (3) can be rewritten as 

nL = (l()[ Fe / H ] - 10[ Fc / H ^) /10F c Wl = 10 [Fc/H]+2.48 _ 10 -0.52^ ^ 

where [Fe/H]p = —3 and [Fe/H]/, = —2.48 are used. The parameters logep, logep, and loge^ 
for the pertinent elements were obtained from the data on two stars with [Fe/H] ps —3 and 
the solar r-process abundances (QW). It was shown that the abundances calculated from the 
above model for a large number of r-process elements were in good agreement with several 
independent data sets over the wide range of —3 < [Fe/H] < — 1. The same agreement was 
also achieved for Sr, Y, Zr, and Ba when the standard solar r-process abundances of these 
elements (Arlandini et al. 1999) were increased by factors of ps 2-6. 
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2.1. Estimates of the P-Inventory and SN 11(H) Yields 

We now extend the above three-component model to the elements below Sr. By the 
assumptions of the model, a fixed amount of ISM dilutes the ejecta from an SN II at [Fe/H] > 
—3. It can be seen from equation (4) that for the standard dilution mass, stars with [Fe/H] ps 
—3 have til ~ 0. So the abundances in these stars only represent mixtures of the P-inventory 
and SN 11(H) contributions: 

10 log6(E) = 10 lo gep (E) +HhX 1Q 1<WE) ( 5 ) 

Thus, with the value of n H obtained from equation (2), the parameters logep(E) and 
\oge H (E) can be determined from the data on any two stars with [Fe/H] ps —3 but dif- 
ferent riH values. 

Figure 2a shows the observed loge values of C to Ge (Westin et al. 2000; McWilliam et 
al. 1995) for HD 122563 ([Fe/H] = -2.74; open squares), HD 115444 ([Fe/H] = -2.99; filled 
circles), and CS 22892-052 ([Fe/H] = -3.03; asterisks) with n H ps 1, 7, and 36, respectively. 
The typical uncertainties in the data shown in Figure 2a are ~ 0.1 dex. It can be seen that 
the abundances of the elements above C are essentially the same for all three stars. In fact, 
when a uniform shift of —0.25 dex is applied to the data on HD 122563 so that all three sets 
of data correspond to essentially the same [Fe/H], significant differences are only found for C 
between the three stars and for Ge between HD 122563 and HD 115444. The differences for 
all the other elements are within the observational uncertainties. As these three stars have 
very different n H values, we conclude that SNe 11(H) produce very little of the elements from 
N to Zn. The abundances of these elements in the three stars are then completely dominated 
by the P-inventory. We take the observed loge values of N to Zn for HD 115444 to be the 
corresponding logep values, thus defining the abundances in the P-inventory. 

An upper limit on log ep for N to Zn can be estimated from 

1 Qloge P (E)+<r(E) > 1 gloge P (E) + ^HD 115444 x 1 Qloge w (E)^ ^ 

where cr(E) ps 0.1 dex is the uncertainty in the observed loge(E) value for HD 115444 and 

n HD115444 _ 7> Equation ( 6 ) gives 

loge H (E)<loge P (E)-1.43 (7) 

for N to Zn. Thus the contributions from an SN II(P) with the standard dilution mass are 
small compared with the existing P-inventory for all the elements between C and Ge. 

The logep values and the upper limits on logep estimated in equation (7) are given 
in Table I. It was argued that the extremely high value of n H ps 36 corresponding to the 
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extremely high r-process abundances observed in CS 22892-052 is due to contamination of 
the surface of this star by the ejecta from the SN 11(H) explosion of a previous massive binary 
companion (Qian & Wasserburg 2001b). The agreement between the observed loge values of 
all the elements from Na to Ni for CS 22892-052 and the corresponding log ep values in Table 
1 further supports the argument that the SN 11(H) yields of these elements are very low. On 
the other hand, the loge(C) value of CS 22892-052 is higher than those of HD 115444 and 
HD 122563 by 1 dex (see Figure 2a). This suggests that a significant amount of C may be 
produced in SNe 11(H) or that C was produced in CS 22892-052 during its evolution. The 
elements C, N, O, Cu, Zn, and Ge will not be discussed further [the evolution of O relative 
to Fe has been discussed in detail in a previous paper (Qian & Wasserburg 2001a)]. Below 
we will focus on the elements from Na to Ni as there are high-quality data on these elements. 

Figure 2b compares the P-inventory (see Table 1) as determined above for [Fe/H] = 
—3 with the observational data at —4 < [Fe/H] < —3. We use the results of Norris, 
Ryan, & Beers (2001) on CD -24°17504 ([Fe/H] = -3.37; squares), CD -38°245 ([Fe/H] = 
-3.98; diamonds), CS 22172-002 ([Fe/H] = -3.61; triangles), CS 22885-096 ([Fe/H] = -3.66; 
asterisks), and CS 22949-037 ([Fe/H] = -3.79; pluses). We have shifted the loge values 
shown for these five stars from the observed values to pass through loge(Fe) = logep(Fe) = 
4.51 so that they are compared at the same Fe abundance. It can be seen that the abundance 
patterns in these stars follow the P-inventory rather closely. Of these stars, CS 22949-037 
exhibits the largest deviation from the P-inventory as it is overabundant in Mg and Si. 
Data on this star from McWilliam et al. (1995) also show that it is overabundant in Na. 
However, the abundance pattern of all the other elements and the ratio Si/Mg for this star 
are essentially the same as for the P-inventory. We thus consider that the P-inventory 
determined by the data at [Fe/H] ps —3 is also an excellent representation of the data at 
—4 < [Fe/H] < —3. We infer that VMS must therefore have a fairly uniform nucleosynthetic 
mechanism if they were responsible for the P-inventory. In §3 we will show that the P- 
inventory must also contain some contributions from SNe 11(H) and ll(L). 



2.2. Estimates of SN II(L) Yields and SN la Yield Pattern 

We now derive the SN ll(L) yields of Na to Ni based on the three-component model. By 
using the logep and loge// values in Table 1 together with equation (1), the loge/, values of 
these elements can be obtained from the observed loge values for HD 178443 with = 2.45 
([Fe/H] = —2.04; McWilliam et al. 1995). The results are given in Table 1. It can be seen 
that typically loge^E) exceeds loge#(E) by ps 2 dex. Thus, although SNe 11(H) are 10 
times more frequent than SNe ll(L) (e.g., QW), numerically the SN 11(H) contributions to 
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the elements from Na to Ni can be ignored in general. 

Using the SN 11(H) and ll(L) yields, we give in Table 1 the fraction a 0jII (E) of the solar 
inventory of an element E contributed by SNe II: 

a 0>n (E) = [nf x 10 lo s^( E ) + nf x io lo s^( E )]/lO lo s e 0( E ), (8) 

where nfj = 10 3 and nf = 10 2 (QW) are the numbers of SNe 11(H) and II(L), respectively, 
that had occurred in a standard dilution mass of ISM prior to solar system formation (SSF). 
The contributions from the P-inventory to the solar inventory are negligible. Table 1 shows 
that SNe II(L) can essentially account for the entire solar inventory of Na, Mg, Si, Ca, Ti, 
and V [a 0) n(E) w 1 for these elements as loge^E) + lognf ~ loge (E)]. By contrast, Table 
1 shows that the a 0) n values for Al, Sc, Cr, Mn, Fe, Co, and Ni only range from 0.13 to 0.47. 
In view of the uncertainties in the data used to derive the loge^ values, we now choose to 
assign the entire solar inventory of Na, Mg, Si, Ca, Ti, and V to SNe II(L) [i.e., we take the 
"corrected" values of logef rr (E) = loge (E) — lognf and ag^E) = 1 for these elements]. 
For Al, Sc, Cr, Mn, Fe, Co, and Ni that have a 0j n values significantly below unity, we 
will use the original calculated values [i.e., loge^ orr (E) = loge^(E) and a jj(E) = a 0) n(E) 
for these elements]. The SN ll(L) yields represented by the loge^ orr values in Table 1 are 
shown in Figure 3 along with the solar abundance pattern that is translated to pass through 
loger r (Mg). 

The majority of the solar inventory of the so-called Fe group elements (Cr, Mn, Fe, Co 
and Ni) may be reasonably attributed to SNe la (e.g. Timmes, Woosley, & Weaver 1995). 
In this case, the part of the solar inventory of an element E in the Fe group contributed by 
SNe la, loge 0) i a (E), can be estimated as 

loge 0)Ia (E) = log[l - a ,„(E)] + loge (E). (9) 

The loge 0i i a values for Cr, Mn, Fe, Co and Ni estimated in equation (9) represent the yield 
pattern of these elements for SNe la and are given in Table I. 

3. The P-Inventory and VMS Yield Pattern 

We have established the P-inventory and SN II(-P) and ll(L) yields of Na to Ni in §2. 
We now focus on the P-inventory. This inventory was considered to be the result of the 
integrated production by VMS over the period prior to the achievement of a metallicity of 
[Fe/H] « —3 in the ISM (WQ). Figure 2b shows that the P-inventory is also an excellent 
representation of the abundance patterns at —4 < [Fe/H] < —3. This suggests that the 
yields of VMS must follow a rather regular pattern. We will seek to establish the yield 
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pattern of VMS from the data at —4 < [Fe/H] < —3. It is important to recognize that 
very old stars with —4 < [Fe/H] < —3 must have low masses of ~ 1 M to have survived 
until the present time. This requires that there be some formation of stars in the normal 
mass range of ~ 1-60 M & in addition to the formation of VMS with masses of > 100 M & 
considerably before a metallicity of [Fe/H] —3 was reached in the ISM. As normal stars 
with masses of ~ 10-60 M are SN II progenitors, we consider that SNe Il(.ff') and H(L) 
must have provided chemical enrichment of the ISM along with VMS at —4 < [Fe/H] < —3. 
Thus, the P-inventory at [Fe/H] = —3 and the abundances in stars with —4 < [Fe/H] < —3 
should be the sum of contributions from VMS and SNe 11(H) and ll(L). The yield pattern 
of VMS must then be obtained by subtracting the SN 11(H) and ll(L) contributions. This 
can be accomplished by using the observed Ba abundances. The observed Ba abundances 
at —4 < [Fe/H] < —3 are particularly important as they also are the key evidence requiring 
a much larger dilution mass in the regime where VMS are active in comparison with the 
regime of [Fe/H] > —3 where VMS are no longer active. 

Unlike SNe II, VMS explode due to pair instability (e.g., Rakavy, Shaviv, & Zinamon 
1967; Bond, Arnett, & Carr 1984; Glatzel, Fricke, & El Eid 1985; Heger & Woosley 2001). 
We do not consider VMS as a source of r-process elements such as Ba. On the other hand, 
small abundances of Ba have been observed at —4 < [Fe/H] < —3 (e.g., Mc William et al. 
1995; McWilliam 1998; Norris et al. 2001). This was previously considered by us to be the 
result of contamination of the VMS-enriched ISM by a small amount of SN 11(H) ejecta that 
contains Ba but very little Fe (WQ). We now consider that some level of "normal" astration 
was present in the regime where VMS were active and that the abundances in the ISM at 
_4 < [Fe/H] < -3 consist of contributions from VMS and SNe 11(H) and II(L). As Ba 
is not considered to be produced by VMS, we may attribute the Ba in a star with these 
metallicities to an effective number hu of SNe II(-P) and an effective number of SNe 
II (L) that had contributed to the ISM from which the star was formed: 

10 lo ge (Ba) = ^ x 10 loge ff (Ba) + ^ x ^^(Ba^ ^ 

where loge^(Ba) = —1.57 or loge L (Ba) = —0.47 corresponds to the abundance of Ba in 
the ISM resulting from a single SN 11(H) or ll(L) for a standard dilution mass of Mf™ « 
3 x 1O 4 M (QW). The lowest observed loge(Ba) value is -2.97 for CS 22885-096 ([Fe/H] = 
—3.66; Norris et al. 2001). Equation (10) gives nn = 0.04 if the Ba abundance in this star 
was provided by a single SN 11(H) or til = 3.16 x 10~ 3 for a single SN II(£). We consider 
that these low values of hu and ul require a much larger dilution mass for the SN II ejecta 
at —4 < [Fe/H] < —3 compared with the standard dilution mass at [Fe/H] > —3. This 
is in accord with the extremely energetic explosions of VMS (e.g., Heger & Woosley 2001) 
that would occur along with SNe II at —4 < [Fe/H] < —3. Specifically, the dilution mass 
M^! MS relevant for these ultra-low metallicities is at least M|^ n /0.04 « 8 x 10 5 M Q and may 
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be as large as M|^ n /(3.16 x 1(T 3 ) f»9x 1O 6 M (the difference in the dilution masses for 
[Fe/H] < —3 and [Fe/H] > —3 was considered earlier by Qian & Wasserburg (2001a) in 
discussing the P-inventory of O). 

We now use the Ba abundances to estimate the SN II contributions for stars with 
—4 < [Fe/H] < —3 in order to obtain the "pure" VMS contributions. As explosions of VMS 
can disrupt regions with a mass much larger than the standard dilution mass for SNe II, 
we assume that the overall mixing occurs over the region governed by VMS explosions. We 
also assume that n H /h L = 10 for a star with —4 < [Fe/H] < —3 and use equation (10) to 
estimate ul as 

n L — ^gloge(Ba)-loge HL (Ba) _ -^glog e(Ba)+0.22 ^ -q 

where we have defined 

IQlogej/iXBa) _ -^q x ^Qloge H (Ba) _|_ ^gloge L (Ba) _ ]^q0.22^ ^2) 

The results presented below are not at all sensitive to the choice of n H /n,L = 10 as the Ba 
yield of SNe II(L) is « 13 times larger than that of SNe 11(H) (see Table 1). For an element 
E other than Ba, the fraction f HL (E) of this element in the star contributed by SNe 11(H) 
and ll(L) is 

f HL (E) =h L x io 1 °s^L(E)-ioge(E) ) (13) 

where log €#l(E) is defined in a similar manner to log e#i(Ba) and can be calculated from the 
loge# and logef rr values in Table 1. The remaining fraction [1 — f HL (E)] is then attributed 
to VMS: 

10 iogevM S (E) = [ X _ f HL (E)] x 10 loge ( E \ (14) 

where loge V Ms(E) represents the total VMS contributions to the element E in the star. The 
above equation also applies to the P-inventory. 

By using the log en and \oge c £ vv values in Table 1 together with equations (11), (13), 
and (14), the fui and logeyMS values for Na to Ni and Sr are calculated from the data on 
four stars with [Fe/H] = —3.98 to —3.61 (Norris et al. 2001) and from the P-inventory as 
determined at [Fe/H] = —3. These results are given in Table 2. The VMS yield patterns 
represented by the log evMS values are also shown in Figure 4. For convenience of comparison, 
the log eyivis values given in Table 2 and shown in Figure 4 have been shifted from the original 
calculated values to pass through the same Fe abundance of logevMs(Fe) = loge^Fe) = 
5.03. It can be seen that the yield patterns of VMS derived from the four stars and that 
derived from the P-inventory are almost identical to each other. The only exceptions are the 
anomalies at Na, Mg, and Si in CS 22949-037 as noted earlier and the varying production of 
Sr. Table 2 also shows that VMS dominated the production of Na to Ni at -4 < [Fe/H] < -3 
and that SNe ll(L) contributed ~ 3-10% of the Fe at these metallicities. 
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The VMS yield patterns derived above may be compared with the SN II (L) yield pattern. 
Figure 4 shows that relative to SNe ll(L), VMS underproduce Na, Al, V, Cr, and Mn, 
overproduce Co, produce a variable amount of Sr, but otherwise have an almost identical 
yield pattern. The Sr production by VMS is small or completely negligible based on the data 
for CD -38°245, CS 22172-002, and CS 22885-096, but may be somewhat larger than that 
by SNe ll(L) based on the data for CS 22949-037 and the P-inventory. Thus there appears 
to be a strong hint for major Sr production by some VMS. 

4. Model for General Chemical Evolution 

We have presented in §§2 and 3 the results on SN 11(H) and ll(L) yields and SN la and 
VMS yield patterns, all of which appear to be rather well-defined. Based on these results, we 
now discuss a model for general chemical evolution of the ISM starting from a zero-metallicity 
state. We will extend the earlier approach for [Fe/H] > —3 (e.g., Qian & Wasserburg 2001a; 
QW) to [Fe/H] < -3. The production by VMS, and SNe 11(H), ll(L), and la, as well as the 
difference in the dilution masses associated with SNe II and VMS, will be taken into account. 
A full approach must address the discreteness of the various production events. We will first 
use a continuous approximation to illustrate some essential features of the model and discuss 
the effects of discrete events later. We use E to represent an element other than Fe. We 
consider the following sequence of evolution for the abundances of E and Fe in the ISM. 
Starting from big bang debris with zero metallicity, i.e., loge(E) = — oo and [Fe/H] = — oo, 
the ISM was first enriched in E and Fe by VMS and SNe 11(H) and ll(L). When a metallicity 
of [Fe/H] « -3 was reached, VMS ceased occurring. At -3 < [Fe/H] < -1, only SNe 11(H) 
and II(L) provided chemical enrichment of the ISM, while at [Fe/H] > -1, SNe 11(H), ll(L), 
and la all contributed. We are concerned here only with the contributions from VMS and 
SNe II and la, but not with the contributions from low mass stars such as asymptotic giant 
branch (AGB) stars. AGB stars would make s-process contributions to the heavy elements 
above Fe at [Fe/H] > — 1. Thus we restrict our discussion to the elements from Na to Ni for 
[Fe/H] < and to the r-process elements (Sr and above) for [Fe/H] < —1. 

4.1. Continuous Model without SN la Contributions 

We here consider a continuous model for evolution of elemental abundances. This model 
is not applicable in regions where the number of nucleosynthetic events is small and addition 
from an individual event causes significant increment in the inventory of an element. This 
approach does provide a good overall view and in particular shows quantitatively that for 
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example, the evolution of Ba relative to Fe with a jump at [Fe/H] « —3 is the result of 
the change from Fe production predominantly by VMS to Fe production by SNe ll(L). The 
continuous model can be discussed without knowing the dilution masses for [Fe/H] < —3 and 
> —3. However, as we have already demonstrated the change in the dilution mass for these 
two regimes in §3, we will incorporate this behavior in the treatment here. The conclusion 
regarding the jump in the evolution of Ba does not depend on the different dilution masses 
for [Fe/H] < -3 and > -3. 

We first discuss the regime of [Fe/H] < —1 that excludes SN la contributions. Consider 
a homogeneous system of gas with a time- dependent total number (H) of H atoms. In the 
continuous approximation, the rate of change in the abundance of E in this gas, d(E/E)/dt, 
is (QW) 

d(E/E) P™S(t) + Pg(t) + P£(t) 

~^T~ ~ (H) ' (i5j 

where P^ MS (t), Pg(t), and P£(t) are the rates for production of E by VMS and SNe 11(H) 
and H(L), respectively. The production rates of individual sources are proportional to the 
absolute yields and the frequencies of these sources. Thus the terms such as P^ MS (t)/(H) 
represent the products of the absolute yields of each event and the frequency of occurrences 
per H atom in the ISM for the individual sources. We consider that the absolute yields of 
VMS and SNe 11(H) and H(L) are fixed. We further assume that the frequencies of these 
sources are proportional to the total number (H) of H atoms in the gas but with different 
proportionality constants for < t < t* and t > t*, where t* is the time when [Fe/H] = —3 
was reached. Then the production rates per H atom, denoted by P^ MS , P^ , and P^ for 
< t < U and by Pg and P^ for t > i*, are all constant. We use a similar notation for Fe. 
Note that Pp e = Pp e = as SNe 11(H) produce no Fe (WQ). The production rates per H 
atom for SNe 11(H) and II (L) at < t < t* (with the tilde symbol) are related to those at 
t > t* (without the tilde symbol) by the very different dilution masses with which the SN II 
ejecta is mixed in these two epochs: 



Pj _ P£ _ Pj _ Mjr 

ri rk Pk m™ s ' 



(16) 



It is found in §3 that the dilution mass for < t < is larger than the dilution mass 

Mf™ 1 for t > U by a factor of ~ 25-300. 

With the notation introduced above, the rates of change in the abundances of E and Fe 
in the gas are 

d(E/H) _ / p™ s + PE + Pi, o<t< u, (17) 



dt \p£ + p£, t> u, 
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and 

d(Fe/H) _ / PJMS + ^ < t < ^ 



I Vg, t > U. ^ 

Eliminating t in equations (17) and (18), we obtain 

gg/H) f (^ MS + ^ + ^)/(^c MS + Pk), [Fe/H] < -3, 

^(Fe/H) \ (Pf + P|)/P£, -3 < [Fe/H] < -1. 1 

Equation (19) with the initial conditions (E/H) = (Fe/H) = at t = can be solved to give 

(F/m = J [( ^ MS + + ^)/(^ e MS + ^ e )](Fe/H), [Fe/H] < -3, 

l/ J 1 (E/H), + [(P# + P|)/^ e ][(Fe/H)-(Fe/H),], -3 < [Fe/H] < -1, 1 J 

where (E/H)* and (Fe/H)* correspond to the abundances of E and Fe at [Fe/H] = —3. 

The result for [Fe/H] < —3 in equation (20) can be rewritten as 

loge(E) = [Fe/H] + loge (Fe) + log[(^ MS + V% + V%) + 7%)]. (21) 

Thus, in the continuous approximation, the evolution of E relative to Fe at [Fe/H] < —3 
follows a straight line of unit slope on a plot of loge(E) as a function of [Fe/H]. The position 
of this line on the plot is determined by the ratio of the relevant total production rates per 
H atom for E and Fe. The result for —3 < [Fe/H] < —1 in equation (20) in the limit of 
(E/H) > (E/H) t and (Fe/H) > (Fe/H) + can be approximated as 



loge(E) » [Fe/H] + log e (Fe)+log[(^ + P|)/P ] 



« [Fe/H] + log e Q (Fe)+\og[(Vg + V£)/Vk\, (22) 

where equation (16) is used in writing the last term. This again represents a straight line 
of unit slope. However, due to the termination of VMS activities at [Fe/H] = —3 and the 
corresponding change in the ratio of the total production rates per H atom for E and Fe, the 
line for equation (22) is offset from that for equation (21) by 



log 



' P, i ' E 1 I Fc + ' Fe 



V VMS +V H + V L I I ^ 



(23) 



Thus the full evolutionary trajectory of loge(E) follows a straight line of unit slope up to 
[Fe/H] = —3 and then asymptotically approaches another straight line of unit slope. The 
sharpness of the transition from the initial line to the asymptotic line can be assessed by the 
slope of the evolutionary trajectory at [Fe/H] = —3: 



dloge(E) 



d[Fe/H] 



d(E/H) 



d(Fe/H) 



(Fe/H)* = / Vi + Vi ) Vl MS + Vk 
(E/H)* [fvMS + fH + f>L \ f>L e 



(24) 
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Note that T-g = for r-process elements such as Ba. For these elements the shift in 
equation (23) is log[(^ c MS + V£ e )/f>£ c ] and the slope in equation (24) is (V^ 1S + Vp c )/Vp e . 
Table 2 shows that SNe ll(L) contributed ~ 3-10% of the Fe in stars with [Fe/H] < —3. So 
(^Ffc MS + ^FeV^Fe ~ 10-30. This results in a rather dramatic jump of ~ 1-1.5 dex between 
the evolutionary tracks of loge(Ba) at [Fe/H] < —3 and [Fe/H] > —2.5. By using the Ba 
and Fe yields of SNe 11(H) and ll(L) in Table 1, the full trajectory of loge(Ba) is calculated 
from equation (20) (see §4.2) and shown in Figure 1 for (V™ s +Vp e ) /Vp e = 10 (solid curve) 
and 20 (long-dashed curve). We recognize that the sparsity of data at [Fe/H] < —3 does 
not permit an accurate determination of the evolution in this region. We consider that 
this evolution may be represented by the corresponding parts of the trajectories shown in 
Figure 1 as the jump between the trends of the data at [Fe/H] < —3 and [Fe/H] > —2.5 is 
adequately reproduced. 



4.2. General Continuous Model 



The above continuous model can be extended to the general case including SN la con- 
tributions to E and Fe in a straightforward manner. For convenience of discussion, we 
designate the following regions for evolution of E relative to Fe: (A) [Fe/H] < —3 with con- 
tributions to E and Fe from VMS and SNe 11(H) and II(L), (B) -3 < [Fe/H] < -2.5 and (C) 
-2.5 < [Fe/H] < -1 with contributions from SNe 11(H) and ll(L), and (D) [Fe/H] > -1 
with contributions from SNe 11(H), ll(L), and la. A typical evolutionary trajectory of 
loge(E) in these four regions is shown in Figure 5a (solid curve). As discussed earlier, the 
trajectory in region A follows a straight line of unit slope determined by the ratio of the total 
rates per H atom for production of E and Fe by VMS and SNe 11(H) and ll(L). Due to the 
termination of VMS activities at [Fe/H] —3 and the corresponding change in the ratio of 
the total production rates per H atom for E and Fe, the trajectory in region B shifts from 
the line in region A. Note that the termination of VMS activities also greatly decreases the 
dilution mass for the SN II ejecta to the standard value. So the loge(E) and [Fe/H] values 
produced by VMS and SNe 11(H) and ll(L) up to [Fe/H] « —3 are soon overwhelmed by 
the additional contributions from SNe 11(H) and II(L) in region B. This typically occurs 
at [Fe/H] ~ —2.5 when a few SNe ll(L) have occurred. The trajectory in region C then 
follows the straight line of unit slope determined by the ratio of the total rates per H atom 
for production of E and Fe by SNe 11(H) and 11(H). At [Fe/H] ~ -1, SNe la began to occur 
with major Fe additions. Thus the ratio of the total production rates per H atom for E and 
Fe may change again at the onset of SN la contributions. The trajectory in region D then 
may shift from the line in region C and reach another asymptotic straight line of unit slope 
after a transition zone. 
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In order to apply the general continuous model to the data on various elements, we need 
to know the ratios of the total production rates per H atom for E and Fe in regions A, B (or 
C), and D: 

7A (E) = (^ MS + ^e +^e)/(^ c MS + ^Fc), (25a) 
7 bc(E) = (Vg + V£)/V£ e = (Vi + V£)/V£ e , (25b) 
7D (E) = (Vg + V^ + V^/iVi + Vk), (25c) 

where V\ and V\ c are the rates per H atom for production of E and Fe, respectively, by 
SNe la. For regions B and C we may estimate 7bc(E) directly from the SN 11(H) and 11(L) 
yields in Table 1 as 

7BC (E) = [10 X 10 1o ^h(E) + 10 lo geL (E)^ 10 lo geL (Fe) ; ^ 

where we have used the frequency ratio of 10:1 for SNe 11(H) and 11(L) (e.g., QW). For 
region A the ratio 7a (E) can be estimated from 

7A (E) _ ( + Vi + V£\f Vi \_ f HL (Fe) 




7bc(E) \ V» + V£ ) \ Vl MS + ) WE) ' 



where ///l(E) and /i?x(Fe) are the fractions of E and Fe at [Fe/H] < —3 contributed by SNe 
11(H) and ll(L) with typical values given in Table 2. 

Now to estimate 7d(E) for region D, we must consider the rates per H atom for Fe 
production by both SNe ll(L) and la, Vp c and Vp e . Once the ratio V\JVy c is determined, 
7 d(E) can be estimated from 

, m tbc(e) + {ri/yhM-Pk/rk) ,„ m 
7D(E) = 1 + (vk/rk) (28) 

by using the SN la yield pattern in Table 1. We recognize that the progenitors of SNe 
11(L) and la evolve very differently and the ratio V\JVy g may not be constant over Galactic 
history. For simplicity, we assume a constant V\JVy c and estimate its value by making use 
of the timescales associated with SN II(L) and la activities. We have assumed that SNe 11(1/ ) 
are responsible for « 1/3 of the solar Fe inventory and SNe la for the remaining fraction of 
w 2/3. Therefore, we have Pp e Ti/(Pp e Tn) 2, where T\ and Tu are the periods of SN la 
and 11(L) production prior to SSF. The ratio T\/T\i is less than unity and depends on the 
time required for a sufficient number of low mass stars to evolve in binaries to provide SN 
la precursors. This time may be ~ 1-3 x 10 9 yr, to be compared with Tu ~ 10 10 yr. We 
consider that 0.7 < T\jT\\ < 1 and use P\J P^ e = 2.5 in equation (28) to estimate 7d(E). 
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Examples of the full evolutionary trajectories for the general continuous model over the 
range of —4 < [Fe/H] < are shown in Figure 6 for Mn and Co along with the available data 
on these two elements. For comparison we also show the same results in the conventional 
[E/Fe] representation. The cusps in the curves at [Fe/H] = —3 and —1 are due to the 
transitions where production in one regime changes to that in the next as described above. 
The intrinsic yield patterns of the different sources are unchanged. Note that the yield 
pattern of Mn, Fe, and Co for VMS is distinct from that for SNe II(L). As VMS turn off and 
SNe II continue with increased production rates per H atom at [Fe/H] = —3, the intrinsic 
differences in the yield patterns of these two sources result in different changes in the ratios 
of the total production rates per H atom and hence, in opposite trends for the evolution 
of Mn and Co in the region of —4 < [Fe/H] < —2.5. These opposite trends as seen in the 
observational data have drawn the attention of Nakamura et al. (1999). They sought to 
explain the observations by changing the yields of the Fe group elements through different 
mass cuts in ab initio models of SN II nucleosynthesis. In our approach, these trends are 
simply explained by the existence of VMS at —4 < [Fe/H] < —3 with a yield pattern of Mn, 
Fe, and Co distinct from that of SNe II(-L). The case of Cr is similar to Mn and can be 
explained in the same manner by our approach. 

The general continuous model can be used to discuss all the other elements in the group 
from Na to Ni as in the above examples for Mn and Co. Summary of the relevant data can 
be found in the figures of McWilliam et al. (1995) and Norris et al. (2001). A common 
feature for all these elements is that the abundance ratios of other elements to Fe become 
rather well defined when [Fe/H] ~ —2.5 is reached. The typical spread in the abundance 
ratio of a particular element to Fe is less than a factor of ~ 3 at a given metallicity in the 
domain of [Fe/H] > -2.5. 

4.3. Effects of Discrete Events 

If there is more than one source for E and Fe as in the case of Ba, the continuous 
approximation is good when a sufficient number of production events of each kind have 
occurred in the ISM. If this requirement is not satisfied, a discrete model should be used. 
The Ba data in Figure 1 show that loge(Ba) increases from « —2.8 at [Fe/H] ~ —4 to ~ —1.8 
at [Fe/H] ~ —3. So the number of Ba producing SNe 11(H) and ll(L) increased by a factor of 
~ 10 over the range of —4 < [Fe/H] < —3 and the continuous approximation should become 
good as [Fe/H] approaches —3. At [Fe/H] > —3, the dilution mass for the SN II ejecta 
decreased greatly to the standard value and a single SN ll(L) would result in [Fe/H] ps —2.5. 
So we expect that the continuous approximation should also hold at [Fe/H] > —2.5 after 
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a few SNe ll(L) had occurred with the standard dilution mass. The evolutionary trends 
exhibited by the Ba data at [Fe/H] < —3 and [Fe/H] > —2.5 are adequately described by 
two straight lines of unit slope offset by ~ 1-1.5 dex as calculated from the continuous model 
(see Figure 1). 

The effects of discrete events are most pronounced in region B of — 3 < [Fe/H] < —2.5 
for r-process elements such as Ba as illustrated schematically in Figure 5b. Due to the 
enormous dilution mass for the SN II ejecta in region A, the typical loge(E) values for 
r-process elements in this region are significantly below the "base level" resulting from a 
single SN 11(H) with the standard dilution mass. This is indicated by the downward arrows 
from the base level in Figure 5b. After the cessation of VMS activities at [Fe/H] « —3, the 
standard dilution mass applies. A single SN 11(H) would then correspond to point a and 10 
SNe 11(H) to point b in region B. As SNe 11(H) produce no Fe, the occurrence of these events 
in region B is represented by the upward arrow from point a to b. A single SN ll(L) with the 
standard dilution mass would correspond to point c at the boundary between regions B and 
C because it produces both Fe and e.g., Ba. As SNe II(-P) are 10 times more frequent than 
SNe ll(L), point d corresponding to 10 SNe 11(H) and a single SN ll(L) would lie on the 
trend line in region C as calculated from the continuous model. For an SN II (L) frequency of 
(10 8 yr) -1 in a standard dilution mass of ISM, region B corresponds to a period of < 10 s yr 
during which SNe 11(H) but not SNe ll(L) would occur with a high probability (WQ). Thus 
the continuous approximation does not apply to region B and the trajectory in this region 
calculated from the continuous model does not provide a physical description of the data on 
e.g., Ba. The large scatter in e.g., loge(Ba) over the narrow range of —3 < [Fe/H] < —2.5 
is accounted for by the occurrence of ~ 1-10 SNe 11(H) that produce Ba but no Fe in the 
discrete model (WQ). 



5. Discussion and Conclusions 

Based on a phenomenological model and the available data on abundances in metal-poor 
stars, we have determined the yields of the elements from Na to Ni for SNe 11(H) and 11(1/ ) 
and the yield patterns of the same elements for SNe la and VMS. In a previous companion 
paper (QW), we determined the yields of SNe 11(H) and ll(L) and the composition of the 
P-inventory for r-process elements (Sr and above). The P-inventory of r-process elements 
is now considered to be the result of concomitant production by SNe 11(H) and H(L) when 
VMS were active. The P-inventory of non-neutron-capture elements (possibly including Sr) 
is considered to be the result of major production by VMS with small contributions from 
SNe 11(H) and H(L). The results presented here and in QW show that the abundances of a 
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large number of elements (of both rapid-neutron-capture and non-neutron- capture origins) 
in any metal-poor star with —3 < [Fe/H] < — 1 can be almost quantitatively determined 
from the observed [Fe/H] and loge(Eu) values for the star by using the three-component 
model including the P-inventory and the contributions from SNe 11(H) and H(L). A similar 
approach also applies to stars with —4 < [Fe/H] < —3, for which the contributions from 
VMS are explicitly included and the contributions from SNe 11(H) and II (L) are determined 
by using the observed loge(Ba) value as the index of r-process production. 

By using the three-component model, it is shown that SNe 11(H) control the abundances 
of all the heavy r-process elements above Ba but contribute almost nothing to the inventory 
of the elements between C and Ge. The solar inventory of all the elements from O to V 
and the solar r-process inventory of all the elements from Sr to Ba are almost quantitatively 
accounted for by SNe H(L). The majority of the solar inventory of the Fe group elements 
(Cr, Mn, Fe, Co, and Ni) is contributed by SNe la. The metal inventory of the universe 
at [Fe/H] < —3 was dominated by the contributions from VMS. Compared with SNe II(L), 
VMS appear to have an almost identical yield pattern but with a distinct deficiency of Na, 
Al, V, Cr, and Mn and a distinct overproduction of Co. All of these discrepant elements 
except for Cr are of odd atomic number. It is found that the metallicity of [Fe/H] « —3 
represents the state when the occurrence of VMS was essentially terminated. While some 
formation of normal stars including SN II progenitors and low mass stars occurred prior to 
the achievement of [Fe/H] « —3, the onset of major formation of normal stars appears to be 
at this metallicity. The rapid rise in the Ba abundance at [Fe/H] —3 (see Figure 1) is due 
to the cessation of VMS activities and the onset of more rapid normal star formation. The 
mass of ISM with which the nucleosynthetic products of VMS and SNe II are mixed in the 
domain of [Fe/H] < —3 is much larger (by a factor of ~ 25-300) than that for SNe II in the 
domain of [Fe/H] > —3 where VMS ceased playing any significant role. 

The above conclusions have implications in four distinct areas: (1) as templates with 
which the results of ab initio nucleosynthetic models of SNe II and VMS may be compared, 
(2) as predictions for the abundances relative to Fe in stars with ultra-low metallicities 
([Fe/H] < —3), (3) as input for discussing the cosmological aspects of condensation and 
dispersion of baryonic matter that led to the formation of galaxies and other aggregations 
of baryonic matter in the universe, and (4) as basis for considering the general chemical 
evolution of the universe. In the following we will present a short summary on each of these 
issues based on the above findings. 
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5.1. Implications for Nucleosynthetic Models of SNe II and VMS 

Based on the three-component model, we have shown in §§2 and 3 that the yield patterns 
of SNe ll{H) and ll(L) and VMS as inferred from the observational data exhibit remarkable 
regularity for a large number of elements. The implications of these regular yield patterns 
for nucleosynthetic models of SNe II and VMS are discussed below. 

5.1.1. SNe 11(H) 

It is commonly thought that in SNe II, the elements below Si are produced by hydro- 
static burning during pre-SN evolution and those above (including Si) are mostly produced 
by explosive burning associated with the shock propagation. In order not to produce a signif- 
icant amount of the elements from Na to Ni, an SN 11(H) must either have a pre-SN structure 
lacking substantial hydrostatic burning shells or for some reason, have all the material below 
the He burning shell fall back onto the central remnant. Stars with masses of ~ 8-11 M 
have very thin shells at the end of their lives and explode due to the collapse of an O-Mg-Ne 
core instead of an Fe core (e.g., Nomoto 1984; Hillebrandt, Nomoto, & Wolff 1984; Mayle & 
Wilson 1988). These stars could be the progenitors for SNe 11(H). Mayle & Wilson (1988) 
calculated that an SN from O-Mg-Ne core collapse would only eject ~ 0.042 M of matter 
and produce ~ 0.002 M of Fe. By assuming that 10 2 SNe II(L) must enrich a standard 
dilution mass of 3 x 10 4 M Q with 1/3 of the solar Fe mass fraction of « 10~ 3 , the 
Fe yield of an SN II(L) is estimated to be « 0.1 M . This is ~ 50 times larger than the 
Fe yield of an SN from O-Mg-Ne core collapse. Given the extremely low total amount of 
ejecta, it is reasonable to expect that SNe from O-Mg-Ne core collapse also do not produce 
any significant amount of other elements below the Fe group. On the other hand, we note 
that neutrino emission from the neutron star produced in any core collapse would drive a 
small amount of material from the neutron star crust. This neutrino-driven wind has been 
suggested as a site of the r-process (e.g., Woosley et al. 1994) although it remains to be seen 
if adequate conditions for the r-process can be obtained in the wind (e.g., Qian & Woosley 
1996; Hoffman, Woosley, & Qian 1997). We consider that SNe from O-Mg-Ne core collapse 
may correspond to SNe 11(H) or at least a subset (see also Wheeler, Cowan, & Hillebrandt 
1998). Further studies of these objects are extremely important as most of the research on 
SNe II tends to focus on progenitors with masses of 12-40 M that lead to explosion from 
Fe core collapse. 

Another possibility to accommodate the nucleosynthetic requirements of SNe 11(H) is 
that SNe II from Fe core collapse would suffer severe fall back of processed material as 
emphasized by Woosley & Weaver (1995). It has been argued that SN 1997D, which ejected 
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only 0.002 M of Fe, was such a case (Benetti et al. 2001). The problem with such a 
scenario for SNe 11(H) is that the severe fall back would be an obstacle for ejection of the 
r-process material from the inner most regions of the SN. In any case, we consider that many 
issues remain to be explored concerning the fall back and possible production and ejection 
of r-process material in SNe from Fe core collapse. 

Given the validity of our model, we require that SNe 11(H) to be the dominant source 
of the heavy r-process elements above Ba but contribute very little of the other elements, 
especially those between C and Ge. We also require SNe II(-H') to be the frequent kind of SNe 
II. Regardless of the theoretical issues involved in addressing the existence of such objects, 
clear predictions can be made for the observational consequences of SNe 11(H). First of all, 
Fe is mostly produced as the radioactive 56 Ni, whose decay powers the light curve of an SN 
II at late stages. The lack of Fe production by SNe 11(H) means that the light curves of such 
events have a very faint tail as in the case of SN 1997D. Based on the frequent occurrences 
of SNe 11(H) in our model, it is likely that future systematic long-term observations of SN II 
light curves may find a large number of cases like SN 1997D. Another observational signature 
of SNe 11(H) relies on the occurrences of these events in binaries. A low mass companion 
would experience surface contamination when an SN II (H) explodes in a binary. As SNe 
11(H) produce mainly the heavy r-process elements above Ba but little else, the surface 
of the low mass companion should be highly enriched in the heavy r-process elements but 
appear more or less normal in e.g., Na to Ni, when compared with other stars with similar 
[Fe/H]. This exactly fits the observations of CS 22892-052 as discussed in §2. We have also 
argued that CS 31082-001 is another example of surface contamination by the SN 11(H) 
ejecta in a binary (Qian & Wasserburg 2001b) based on the highly enriched abundances of 
heavy r-process elements but very low value of [Fe/H] = —2.9 for this star (Cayrel et al. 
2001). There are no available data on the elements in the group of Na to Ni other than 
Fe for CS 31082-001. As the value of [Fe/H] = -2.9 is almost the same as that for the P- 
inventory, we predict that future observations of this star should detect abundances similar 
to the P-inventory for all the other elements in the group of Na to Ni. The abundances of 
heavy r-process elements observed in CS 31082-001 are extremely high and correspond to 
n H ~ 100. Thus the observations of Na to Ni in this star may be used to set a severe limit on 
the production of these elements by SNe 11(H). We urge that such observations be carried 
out soon. 
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5.1.2. SNell(L) 

The SN 11(1/ ) yields presented here are inferred from the data on a star with [Fe/H] = 
—2.04 corresponding to several contributing SNe II(L). These yields almost exactly follow 
the solar abundance pattern up through V. This suggests that the SN 11(L) yields may be 
interpreted as the contributions from a single type of event and that the yield pattern is 
constant, i.e., independent of metallicity over Galactic history. The extremely regular SN 
II (L) yields presented here may be compared with the yields calculated from nucleosynthetic 
models of SNe II, for which diversity appears to be the rule. Figure 7 compares the SN ll(L) 
yields (filled circles connected by solid curves) with the yield patterns calculated from models 
(Woosley & Weaver 1995) of SNe II with 20 M progenitors of solar (triangles) and 0.01 times 
solar (squares) metallicities. The model yields have been shifted to pass through the SN II (L) 
yield of Mg in panel a and that of Fe in panel b. These two choices of the reference point are 
made to reflect that the elements below and above Si have different nucleosynthetic origins 
in SNe II. It can be seen that the model yields are dependent on metallicity and that both 
model yield patterns only resemble the SN ll(L) yield pattern in local regions. 

To emphasize the regularity of SN II (L) yields inferred from observations, we show in 
Figure 8 the data on log(Si/Mg) over the wide range of —4 < [Fe/H] < —1.3. It can be seen 
that most 75%) of the data are within a factor of 2 of Si/Mg = 1 ^ (Si/Mg) and that 
there are no evidence for high Si/Mg ratios of > 3. In contrast to this regularity, the values 
of Si/Mg in most models (e.g., Woosley & Weaver 1995) vary over a wide range of ~ 0.1-10. 
This variability is expected as Mg is produced by hydrostatic burning but Si by a mixture of 
hydrostatic and explosive burning. We emphasize that averaging over many SNe II does not 
provide a viable way to reconcile the diverse Si/Mg ratios in the models with the uniform 
Si/Mg values exhibited by the data. This is because the same uniformity of Si/Mg applies 
to the data at [Fe/H] ~ —2.5, which corresponds to Fe contributions from few SNe II(-L) . 

The modelling of explosive burning is extremely sensitive to many uncertainties in our 
understanding of SNe II as emphasized by Woosley & Weaver (1995) and Thielemann, 
Nomoto, & Hashimoto (1996). This sensitivity is commonly explored by adopting vari- 
ous mass cuts for the SN II ejecta. However, introduction of mass cuts deviates from the ab 
initio approach and must rely on empirical guidance such as the Fe yields inferred from SN 
II light curves. Note that the typical inferred Fe yields are ~ 0.1 M (e.g., Thielemann et 
al. 1996), which are consistent with the Fe yield estimated for SNe ll(L) in the preceding 
discussion on SNe II(^) . It is conceivable that by choosing a mass cut to give the empirical 
Fe yield, SN II models would give more convergent yield patterns similar to the SN II (L) 
yield pattern derived here from the observed abundances in metal-poor stars. A deeper issue 
then arises concerning the physics and astrophysics underlying the mass cut. 
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5.1.3. VMS 

We have shown that the yield pattern of VMS is also quite regular. In particular, this 
pattern is almost identical to the SN II (L) yield pattern except for the distinct deficiency 
at Na, Al, V, Cr, and Mn, the overproduction of Co, and the variable production of Sr (see 
Figure 4). We will not discuss Sr any further. All the other discrepant elements are of odd 
atomic numbers except for Cr. The underproduction of the elements of odd atomic numbers 
by VMS is expected from the requirement of neutron excess to produce these elements. 
During the evolution of an SN II progenitor, significant weak interaction such as [3 + decay 
and electron capture would take place during hydrostatic burning to increase the neutron 
excess. However, VMS encounter pair-instability after core He burning and there is no stage 
of stable post-He burning to increase the neutron excess (Heger & Woosley 2001). Thus, 
VMS are expected to greatly underproduce the elements of odd atomic numbers relative to 
SNe II. However, the underproduction of Cr of even atomic number, the overproduction of 
Co of odd atomic number, and the similar production of Sc of odd atomic number by VMS 
relative to SNe II (L) as found here are not readily explained. 

Figure 9 compares the VMS yield patterns inferred from the data on four stars with 
—4 < [Fe/H] < —3 and that derived from the P-inventory with the yields calculated from 
two limiting models (Heger & Woosley 2001) of VMS with He core masses of 70 M & (dashed 
curve) and 130 M (solid curve). The model yields have been shifted to pass through 
loge(Mg) = 5.58, the location of most logeyMs(Mg) values, in panel a and loge(Fe) = 5.03, 
the location of all logevMs(Fe) values, in panel b. These two choices of the reference point 
are made again to reflect the relative robustness of the yields for Mg and Fe. It can be seen 
that below Ti, the model yields are quite robust and are a fair representation of the elements 
of even atomic numbers. By contrast, the production of the elements above Ti in the models 
increases with the He core mass. It is also evident that the production of the elements of 
odd atomic numbers and Cr cannot be simply accounted for by the models in a quantitative 
manner. Comparison with models for intermediate He core masses does not change the basic 
conclusions. 

The explosion energy in the VMS models of Heger & Woosley (2001) lies in the range of 
~ 10 52 -10 53 erg and increases with the VMS mass. These values are larger than the typical 
SN II explosion energy of « 10 51 erg by factors of ~ 10-100. From detailed numerical studies, 
Thornton et al. (1998) found that the total mass swept up by an SN remnant scales with the 
explosion energy to a power of 6/7. This suggests a dilution mass for VMS that is ~ 7-50 
times the standard dilution mass of M^ u « 3 x 10 4 M Q for SNe II and is consistent with the 
lower end of the enormous dilution masses of ~ 1O 6 -1O 7 M derived here for VMS. 

Assuming that a single VMS gives [Fe/H] ~ — 4 in a dilution mass of ~ 10 6 M , we estimate 
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that the Fe yield of VMS is ~ 0.1 M . This also implies that ~ 10 VMS had occurred in this 
dilution mass prior to the cessation of VMS activities at [Fe/H] ~ —3 [Qian & Wasserburg 
(2001a) concluded that many VMS are required to provide the P-inventory of O]. 

5.1.4. Summary 

From the detailed nucleosynthetic models of SNe II and VMS we would expect there 
to be a wide range in the relative production of the elements from Na to Ni. However, 
we have shown in §4 that the evolution of the other elements relative to Fe over a wide 
range in [Fe/H] can be adequately described by assuming fixed yield patterns of SNe II (L) 
and VMS (see Figure 6). This poses a serious dilemma. One obvious solution would be 
that the sampling of the ISM by each star represents contributions from many SNe II and 
VMS whose average is strongly convergent. We well recognize that mixing can hide grossly 
diverse sources. However, it is not likely that the regular behavior exhibited by the data 
over the wide range in [Fe/H] can be all attributed to mixing. For example, in the domain of 
—2.5 < [Fe/H] < —2, the abundances of Na to Ni should be dominated by the contributions 
from only several SNe I1(L) and mixing would not be effective. But the data in Figure 
6 show that the relative abundances of Mn and Co to Fe already are well defined in this 
domain. The data on other elements show the same behavior and can be found in the figures 
of McWilliam et al. (1995) and Norris et al. (2001). 

While a statistical equilibrium among the a-nuclei from 28 Si to e.g., 52 Fe, can be possibly 
at work to set a fixed yield pattern of the corresponding elements, there is no explanation, 
to our knowledge, of the yield patterns outlined here for SNe 11(1/ ) and VMS. The observed 
constancy of Si/Mg « 1 over a wide range in [Fe/H] is in sharp contrast to what is found 
in models of SNe II and VMS. We are pressed to inquire whether or not a new paradigm 
for SN II and VMS models should be explored that would treat the possibility of much 
more restrictive evolutionary paths for SNe II and VMS. From the arguments presented 
here, the major efforts to model the solar abundance pattern by many sources with widely- 
varying yield patterns have not led to a satisfactory conclusion. We recognize that ab initio 
calculations are both complex and difficult. It is possible that some unspecified dynamical 
effects involving partially quenched evolution through high temperature and high density 
stages may be of importance in the production of Mg and Si. The phenomenological approach 
used here can in no way describe the detailed nuclear astrophysical evolution of SNe II or 
VMS. However, we consider that the rather robust yield patterns of SNe II(L) and VMS 
are quite directly obtained from the observational data and we cannot identify any obvious 
errors that would change the results. The only substantial variation is found for the VMS 
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yields of Na, Mg, and Si from the data on one star (CS 22949-037). This lack of variation in 
SN ll(L) and VMS yields is remarkable in view of the different data sources used. We do not 
consider that this is plausibly attributable to some grand averaging. Instead, there may be 
possible mechanisms that would give sharply convergent yield patterns nearly independent 
of e.g., the mass or metallicity of SN II progenitors. Certainly the observational data should 
serve as a guide to theoretical exploration of such mechanisms. 



5.2. Implications for Abundances at Ultra-low Metallicities 

Based on the model presented here, the abundances at ultra-low metallicities ([Fe/H] < 
—3) can be described in terms of the yield pattern of VMS and the yields of SNe 11(H) 
and ll(L). For the elements from Na to Ni, the dominant contributions are from VMS. 
Thus, to first approximation, the observed abundances at [Fe/H] < —3 should follow the 
VMS yield pattern. A more detailed study can be carried out by using the Ba abundance 
as the index for r-process production to subtract the concomitant contributions from SNe 
11(H) and ll(L). Such a study may reveal more clearly the features at the elements of 
odd atomic numbers and Cr as these elements receive more contributions from SNe H(L). 
All the r-process elements above Sr should be attributed to SNe 11(H) and ll(L). The 
abundances of these elements at [Fe/H] < —3 may be complicated by the effects of discrete 
events and issues of the dilution mass. For example, a low mass star might form in a 
region where an SN H(-f^) just occurred but no nearby VMS explosion occurred yet. In this 
case, the star would sample the enrichments by an SN 11(H) with a standard dilution mass 
and thus show significant enhancements in the heavy r-process elements while maintaining 
an ultra-low [Fe/H]. Note that a star sampling the enrichments by an SN II(L) with a 
standard dilution mass would move out of the domain of ultra-low metallicities as it would 
have [Fe/H] R3 —2.5 (see Figure 5b). As SNe 11(H) are frequent, future observations at 
[Fe/H] < —3 may find stars with unusually high abundances of heavy r-process elements. 
In any case, given sufficient sensitivity, future observations may detect the heavy r-process 
elements such as Eu that were produced by SNe 11(H) concurrent with VMS activities at 
ultra-low metallicities. It is evident that there is an urgent need for more extensive and 
precise observations of abundances in stars with metallicities well below —3. 
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5.3. Implications for Formation of Galaxies and Chemical Evolution of the 

Universe 

In considering broader implications of the results presented here, it should be noted 
that the conclusions are only based on observational data and a phenomenological model. 
No direct considerations are made of cosmological models or the dynamics and mechanisms 
associated with SN and VMS evolution. Time is, for the most part, only considered in terms 
of sequence. Increments in the inventory of elements associated with newer generations of 
stars are taken as the measure of time sequence. The other guiding rules are the general 
aspects of stellar evolution. 

The following general evolution sequence is inferred: (1) formation of VMS made of 
big bang debris, (2) concomitant formation at some level of a normal stellar population, 
(3) "metal" enrichment with overwhelming contributions from VMS until a metallicity of 
[Fe/H] —3 was reached, (4) cessation of VMS activities and onset of major formation of 
normal stellar populations as marked by a sharp increase of contributions from the frequent 
SNe 11(H) at [Fe/H] ps —3, (5) continued chemical enrichment with dominant contributions 
from SNe 11(H) and ll(L) until [Fe/H] ~ —1 was reached, (6) onset of SN la contributions 
at [Fe/H] ~ — 1 when a sufficient number of low mass stars had evolved in binaries to serve 
as precursors for SNe la, and (7) continued chemical enrichment with contributions from 
SNe 11(H), ll(L), la and other evolved low mass stars. In the regime of [Fe/H] < —3 where 
VMS play a predominant role, we infer that in aggregations of matter where astration has 
begun, the explosion of a VMS redistributes and mixes matter over a mass scale of ~ 10 6 - 
1O 7 M . This involves fresh nucleosynthetic products of the VMS and any SN II debris in 
the region. Thus, when VMS activities cease at [Fe/H] ps —3, the dilution mass with which 
nucleosynthetic products are mixed drops by a few orders of mag nitude from ~ 1O 6 -1O 7 M 
for [Fe/H] < -3 to the value of ps 3 x 10 4 M associated with SNe II for [Fe/H] > -3. 

In typical cosmological models of hierarchical structure formation, the condensation of 
dark matter is considered to provide the gravitational potential wells in which some baryonic 
matter collects. The typical binding energies of these wells correspond to escape velocities 
of ~ 400 km s -1 . The VMS explosions are more than sufficient to disrupt the baryonic 
condensates in such wells. We consider that in the regime of [Fe/H] < —3, the increases in 
"metallicity" are due to the sequence of local condensation of baryonic matter, formation 
of VMS and some normal stars, and dispersion by VMS explosion. There are then later 
reaggregations of some baryonic matter and further enrichments through repeating the above 
sequence. The cessation of VMS activities is due to the increase in metallicity that favors 
formation of stars with a more normal distribution in mass. This permits the aggregation of 
some matter without disruption. Subsequent to the cessation of VMS activities at [Fe/H] ps 
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—3, large scale dispersion of baryonic matter becomes very limited and large aggregation 
of baryonic matter that leads to normal galactic structures and galactic chemical evolution 
begins to occur. A recent study by Bromm et al. (2001) indicates that substantial formation 
of a normal stellar population could only occur when the metallicity is substantially above 

5 x 10 -4 times the solar value. 

In considering the domain where VMS and SNe II coexist (region A in Figure 5), it is 
necessary to discuss the timing of events. Based on the difficulty of forming stars with lower 
masses from big bang debris (e.g., Bromm et al. 2001), the first objects formed would have 
to be VMS. Upon explosion these would disrupt the original baryonic aggregates formed 
in the potential wells of dark matter. Reaggregation and repetition of such events would 
increase the metallicity to allow formation of both VMS and lower mass stars. Certainly 
some lower mass stars began to form by [Fe/H] ~ —4. Consider a "cloud" which contains 
VMS and SN II progenitors. The VMS would explode first and disrupt the cloud. The 
subsequent SN II explosions would disperse ejecta into a hot and tenuous medium that does 
not permit any immediate star formation. After sufficient cooling, the ejecta from the VMS 
and SNe II and the general medium would reaggregate under the influence of the dark matter 
potential. During this process, the nucleosynthetic products of VMS and SNe II would be 
effectively mixed over the entire mass of the baryonic aggregate prior to astration. Thus, 
the effective dilution mass for SNe II (much larger than the standard dilution mass) is the 
same as that for VMS in the regime where VMS are active. When [Fe/H] ks —3 is reached, 
astration into lower mass stars becomes very efficient and formation of VMS is truncated. 
The total timescale over which normal astration becomes dominant is not known. Based on 
considerations of data on damped Lja systems, Wasserburg & Qian (2000b) estimated that 
this timescale was typically several 10 9 yr after the big bang. The question of reaggregation 
is a complex one as it involves condensation from ionized matter. It is known that most of 
the baryonic matter resides in the ionized intergalactic medium (IGM). This matter has not 
yet formed and may never form galaxies. 

Observations of damped Lja systems over a wide range in redshift (z ~ 1.5-4.5) show 
that the lowest [Fe/H] observed is ~ —2.7 (e.g., Prochaska & Wolfe 2000; Prochaska, Gawiser, 

6 Wolfe 2001). The damped Lja systems are thought to represent proto-galaxies. The lower 
bound on [Fe/H] in these systems was interpreted to reflect the transition to normal astration 
at [Fe/H] « —3 over an extended timescale after the big bang (Wasserburg & Qian 2000b). 
As this is the same effective bound for a wide range of z, it follows that subsequent to 
the cessation of VMS activities at [Fe/H] —3, most of the baryonic matter remained 
dispersed in the universe to serve as a reservoir for formation of proto-galaxies. Based on 
the model presented here, this matter should exhibit the chemical enrichments corresponding 
to the P-inventory that resulted from the integrated production by VMS and some SNe II 
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prior to the achievement of [Fe/H] « —3. Assuming the yields of Heger & Woosley (2001), 
Oh et al. (2001) have shown that a fraction ~ 10 -5 -10~ 4 of all baryonic matter must be 
processed through VMS in order to account for the P-inventory at [Fe/H] —3. Using 
the VMS luminosities of Bromm, Kudritzki, & Loeb (2001), they have further argued that 
this amount of processing would provide ~ 10 photons per baryon at energies sufficient to 
ionize H and He. This suggests that VMS may be sufficient to explain the Gunn-Peterson 
effect (Gunn & Peterson 1965), which requires that most of the baryonic matter be ionized. 
How effective the ionization by VMS may be for [Fe/H] substantially below —3 is not yet 
explored. The P-inventory given here is thus considered to represent the abundances in 
dispersed ionized baryonic matter and should be compared with what is observed in e.g., the 
1GM. If an acceleration mechanism were to exist in the dispersed ionized baryonic matter, a 
cosmic ray component might be produced with a composition reflecting the P-inventory. 

It has long been recognized that the ratio E/Fe for the a-elements such as Mg, Si, and 
Ca in the early Galaxy is higher than the solar value by a factor of ~ 3. This was due to the 
additional Fe production by SNe la at later times (e.g., Tinsley 1980). Only 1/3 of the 
solar Fe inventory was produced by SNe II. The ratio E/Fe for the a-elements in the early 
Galaxy was considered to reflect the production of these elements by SNe II. We have shown 
that the production of these elements at [Fe/H] < —3 was dominated by VMS. However, the 
yield pattern of Mg, Si, Ca, Ti, and Fe for VMS is almost identical to that for SNe H(L). 
Thus, the ratio E/Fe for the ct-elements stays approximately constant prior to the onset of 
SN la contributions although the dominant production sources have changed at [Fe/H] —3 
due to the cessation of VMS activities. The special status of [Fe/H] « —3 is demonstrated by 
the sharp increase in Ba abundance at this metallicity (the same sharp increase also occurs 
at e.g., [Si/H] « -2.5 as [Si/Fe] « 0.5 at low [Fe/H]; see Figure 2 of Oh et al. 2001). The 
cosmological epoch prior to the achievement of [Fe/H] —3 remains to be explored. 
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Fig. I. — Data (filled circles: Ryan, Norris, & Beers 1996; Norris et al. 2001; filled diamonds: 
McWilliam et al. 1995; McWilliam 1998; open squares: Burris et al. 2000; lower and upper 
open triangles: HD 122563 and HD 115444, Westin et al. 2000; open circle: CS 22892-052, 
Sneden et al. 2000) on loge(Ba) as a function of [Fe/H]. The Ba yield of a single SN II(ff) 
is shown as the dotted line and that of a single SN 11(f) is 1.1 dex above this. The straight 
dot-dashed line represents the evolutionary trend for production only by SNe II (if) and II (L) 
with a frequency ratio of 10:1. The extension of this line below [Fe/H] ~ —2.5 corresponds 
to greater than standard dilution masses but no other Ba or Fe sources. The solid and 
long-dashed curves are for the continuous model of evolution with VMS activities and the 
associated large dilution mass at [Fe/H] < —3. It is considered that VMS produce no Ba 
but dominated the Fe production until the cessation of VMS activities at [Fe/H] —3. The 
solid and long-dashed curves assume that SNe 11(f) contributed 10% and 5%, respectively, 
of the Fe at [Fe/H] < —3. The continuous approximation fails and a discrete model must be 
used in the transition region between the two short-dashed lines (see discussion of Figure 5). 

Fig. 2.— (a) Data on three stars (open squares: HD 122563 with [Fe/H] = -2.74, filled 
circles connected with solid curves: HD 115444 with [Fe/H] = —2.99, Westin et al. 2000; 
asterisks: CS 22892-052 with [Fe/H] = -3.03, McWilliam et al. 1995) that are used to infer 
the P-inventory and SN II (if) yields. Note the remarkable similarity in the abundances of all 
the elements above C for these stars with widely varying Eu abundances corresponding to SN 
Il(ff) contributions for nn ~ 1, 7 and 36. (b) The P-inventory (filled circles connected with 
solid curves) calculated from the data in (a). The data on five stars (squares: CD — 24°17504 
with [Fe/H] = -3.37 and no Sr data, diamonds: CD -38°245 with [Fe/H] = -3.98, triangles: 
CS 22172-002, asterisks: CS 22885-096, pluses: CS 22949-037; McWilliam et al. 1995 for 
Na and Norris et al. 2001 for all the other elements) are shown for comparison. These data 
have been shifted from the observed abundances to pass through the same Fe abundance of 
loge(Fe) = logep(Fe) = 4.51 and are almost indistinguishable from the P-inventory. There 
appear to be some anomalies at Na, Mg, and Si for one star (pluses: CS 22949-037), some 
spread in Mn, and a wide scatter in Sr. 

Fig. 3. — The SN 11(f) yields (filled circles connected with solid curves) calculated from the 
data on HD 178443 (McWilliam et al. 1995) by subtracting the small contributions from 
the P-inventory and SNe II(ff). The solar abundance pattern translated to pass through 
loge(Mg) = loge^(Mg) = 5.58 (dashed curve) is shown for comparison. With the exception 
of Al and Sc, the translated solar abundance pattern up through V coincides with the 
SN 11(f) yields. The differences for the Fe group elements (Cr, Mn, Fe, Co, and Ni) are 
attributed to the sum of SN la contributions corresponding to the loge 0) i a values in Table 1. 

Fig. 4. — The VMS yield patterns calculated from the data on four stars with [Fe/H] < —3 
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(diamonds: CD -38°245, triangles: CS 22172-002, asterisks: CS 22885-096, pluses: CS 
22949-037; McWilliam et al. 1995 for Na and Norris et al. 2001 for all the other elements) 
and that calculated from the P- inventory at [Fe/H] = —3 (open circles) by using Ba as 
the index of SN II production. These patterns have been translated to pass through the 
same Fe abundance of loge(Fe) = loge^Fe) = 5.03 for comparison with the SN II(X) yield 
pattern (filled circles connected by solid curves). Note the general similarity between the 
VMS and SN ll(L) yield patterns. However, VMS show some high production of Na, Mg, 
and Si from the data on one star (pluses: CS 22949-037), and in all cases a significantly lower 
production of Na, Al, V, Cr, and Mn and a significantly higher production of Co. Except 
for CS 22885-096 (asterisks with no Sr assignment), the Sr in all the other three stars and 
in the P-inventory is not explained by SN II coproduction with Ba. This is considered as a 
strong hint of Sr production by some VMS. 

Fig. 5. — (a) Schematic illustration of continuous evolution of loge(E) as a function of [Fe/H] 
in four different regions (A, B, C, and D). The production sources in each region are indicated. 
The cusps at [Fe/H] —3 and ~ —1 are due to the transitions where production in one region 
changes to that in the next (see text), (b) Schematic illustration of the effects of discrete 
events on the evolution of r-process elements such as Ba. Due to the enormous dilution mass 
associated with VMS activities in region A, the typical abundances in this region are well 
below the base level resulting from a single SN ll{H) with the standard dilution mass. The 
VMS activities cease at [Fe/H] « —3 and the dilution mass drops sharply to the standard 
value. Point a in region B corresponds to a single SN 11(H) and point b to 10 SNe 11(H). As 
SNe 11(H) produce no Fe, the occurrence of these events in region B follow the upward arrow 
from point a to b. Point c corresponds to a single SN II(L) and point d to 10 SNe 11(H) and 
a single SN ll(L), all with the standard dilution mass. The continuous approximation fails 
and a discrete model must be used in region B (see text). 

Fig. 6. — (a) Evolution of Mn relative to Fe for the continuous model (solid curves) is 
compared with the data (filled circles: Ryan et al. 1996; Norris et al. 2001; filled diamonds: 
McWilliam et al. 1995; open squares: Gratton 1989; open triangles: Westin et al. 2000). The 
asterisk corresponds to the solar data. The parameters used are /#z,(Mn) = 0.3, /ffi(Fe) = 
0.1, and Vp c /Vp c = 2.5. No attempt is made to fit the model to the data, (b) Same as (a) 
but for Co (with open squares indicating data from Gratton & Sneden 1988 and with only 
one different parameter //jx(Co) = 0.01 for the model). 

Fig. 7. — The SN ll(L) yield pattern (filled circles connected by solid curves) is compared 
with the yield patterns calculated from models (Woosley & Weaver 1995) of SNe II with 
20 M & progenitors of solar (triangles) and 0.01 times solar (squares) metallicities. The model 
yields have been shifted to pass through loge(Mg) = logez,(Mg) = 5.58 in panel a and 
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loge(Fe) = loge_L(Fe) = 5.03 in panel b. The solar abundance pattern translated to pass 
through loge(Mg) = 5.58 (dashed curve) is also shown in panel a. 

Fig. 8. — The observed log(Si/Mg) values over the range of —4 < [Fe/H] < —1.3 (filled 
circles: Ryan et al. 1996; Norris et al. 2001; filled diamonds: McWilliam et al. 1995; 
open squares: Gratton & Sneden 1988; open triangles: Westin et al. 2000). The solid line 
represents the solar value of log (Si/Mg) Q = —0.03. Note that almost all (rs 75%) of the 
data lie within ±0.3 dex (between the dashed lines) of log (Si/Mg) = 0. There is no evidence 
of very high Si/Mg ratios. 

Fig. 9. — The VMS yield patterns derived from the data on four stars with [Fe/H] < —3 
and that derived from the P-inventory at [Fe/H] = —3 (same symbols as in Figure 4) are 
compared with the yield patterns calculated from models (Heger & Woosley 2001) of VMS 
with He core masses of Mn e = 70 M & (dashed curve) and 130 M & (solid curve). The model 
yields have been shifted to pass through loge(Mg) = 5.58 in panel a and loge(Fe) = 5.03 in 
panel b. 
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Table 1. P-Inventory, SN 11(H) and II (L) Yields, and SN la Yield Pattern 





E 


logep(E) 


loge H (E) 


log£i(E) 


loge (E) 


"0,ii (E) 


loge«> rr (E) 


°0,n( E ) 


loge ,i a (E) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Group 


1: 


















6 




< 5.46 






8.56 








— oo 


7 


N 


6.30 


< 4.87 




8.05 








— oo 


8 


u 


6.60 


~ 5.17 


6.93 


8.93 


1.00 


6.93 


1.00 


— oo 


Group 


2: 


















11 


Na 


3.34 


~ 1-91 


4.24 


6.33 


0.81 


4.33 


1.00 


— oc 


12 


Mg 


5.13 


~ 3.70 


5.62 


7.58 


1.09 


5.58 


1.00 


— oo 


13 


Al 


3.12 


~ 1-69 


3.90 


6.47 


0.27 


3.90 


0.27 


— oo 


14 


Si 


5.02 


<3.59 


5.47 


7.55 


0.84 


5.55 


1.00 


— oo 


20 


Ca 


3.75 


< 2.32 


4.40 


6.36 


1.10 


4.36 


1.00 


— oo 


21 


Sc 


0.28 


< -1.15 


0.77 


3.10 


0.46 


0.77 


0.46 


— oo 


22 


Ti 


2.43 


< 1.00 


2.87 


4.99 


0.76 


2.99 


1.00 


— oo 


23 


V 


1.15 


< -0.28 


1.87 


4.00 


0.74 


2.00 


1.00 


— oc 


24 


Cr 


2.40 


< 0.97 


3.17 


5.67 


0.32 


3.17 


0.32 


5.50 


25 


Mn 


1.90 


<0.47 


2.56 


5.39 


0.15 


2.56 


0.15 


5.32 


26 


Fe 


4.51 


< 3.08 


5.03 


7.51 


0.33 


5.03 


0.33 


7.34 


27 


Co 


2.24 


< 0.81 


2.05 


4.92 


0.13 


2.05 


0.13 


4.86 


28 


Ni 


3.25 


< 1.82 


3.50 


6.25 


0.18 


3.50 


0.18 


6.16 


Group 


3: 


















29 


Cu 


0.53 


< -0.90 


? 


4.21 


? 


? 


? 


? 


30 


Zn 


1.86 


< 0.43 


? 


4.60 


? 


? 


? 


? 


32 


Ge 


< -0.05 


? 


? 


3.41 


? 


? 


? 


— oo 


Group 


4: 


















38 


Sr 


0.13 


-1.30 


0.35 


2.90 


0.34 


0.35 


0.34 


— oo 


56 


Ba 


ps -1.80 


-1.57 


-0.47 


2.13 


0.45 


-0.47 


0.45 


— oc 



Note. — Results are given for four groups of elements with the focus on groups 2 and 4. The P-inventory in col. 
(3) and SN Il(-ff) yields in col. (4) are calculated from the data shown in Figure 2a for groups 1-3. The SN II(L) 
yields in col. (5) are calculated from col. (3) and the data on HD 178443 (McWilliam ct al. 1995) for group 2 and 
that for O is calculated by attributing the solar O inventory to 10 2 SNe II(L). Col. (6) gives the solar inventory 
as measured in the photosphere (Anders & Grevesse 1989). Col. (7) gives the fraction of the solar inventory of an 
clement E contributed by SNe II as calculated from cols. (4)-(6). The corrected values in cols. (8) and (9) only 
reflect small adjustments for those elements with aQ ii(E) 1. The SN la yield pattern is represented by the part 
of the solar inventory contributed by SNe la in column (10) as calculated from cols. (6) and (9). The results for 
group 4 are taken from QW but with logep(Ba) estimated from the data at [Fe/H] < —3 (see Figure 1). 

a Atomic number. 



Table 2. VMS Yield Patterns 
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CD 


-38°245 


CS 22172-002 


CS 22885-096 


CS 22949-037 


P-inventory 


/fli(E) 


log evms(E) 


JHL (&) 


1/~*<T Ci i» »f-i I TV, 1 

lug &VMS K 1 ^) 




lug fcVMSV 1 ^/ 


JHL(£J) 


log fcvMS K 1 ^ ) 


jhl y^) 


10 6 fc VMS V 11 '/ 


"1 "1 


in a 




o.Oo 






0.14 


3.56 


< 0.01 


5.87 


0.26 


3.77 


12 


Mg 


0.11 


5.55 


0.07 


5.46 


0.02 


5.62 


0.02 


6.36 


0.07 


5.66 


13 


Al 


0.49 


2.97 


0.36 


2.88 


0.13 


3.16 


0.23 


3.49 


0.16 


3.60 


14 


Si 


0.15 


5.34 


0.06 


5.44 


0.03 


5.51 


0.03 


6.15 


0.09 


5.54 


20 


Ca 


0.22 
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0.11 


3.99 


0.04 


4.15 


0.11 


4.33 


0.11 


4.26 


21 


Sc 


0.12 


0.69 


0.07 


0.64 


0.03 


0.77 


0.11 


0.75 


0.08 


0.80 


22 


Ti 


0.14 


2.82 


0.08 


2.80 


0.05 


2.67 


0.14 


2.84 


0.10 


2.95 


23 
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> 0.09 


< 2.04 






> 0.04 


< 1.82 


> 0.07 


< 2.19 


0.19 


1.62 


24 


Cr 


0.44 
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0.21 


2.48 


0.10 


2.57 


0.36 


2.49 


0.15 


2.89 


25 


Mn 


0.51 


1.59 


0.29 


1.69 


0.07 


2.10 


0.40 


1.82 


0.12 


2.40 


26 


Fe 


0.10 


5.03 


0.05 


5.03 


0.03 


5.03 


0.11 


5.03 


0.09 


5.03 


27 


Co 


0.01 


2.97 


0.01 


3.03 


< 0.01 


3.04 


0.01 


3.06 


0.02 


2.79 


28 


Ni 


0.07 


3.66 


0.03 


3.80 


0.01 


3.96 


0.05 


3.86 


0.05 


3.79 


38 


Sr 


0.54 


-0.57 


0.85 


-1.58 


1.10 


— oo 


0.09 


0.53 


0.05 


0.67 



Note. — Results are calculated from the data on the four listed stars with [Fe/H] < —3 (McWilliam et al. 1995 for Na; Norris et al. 2001 for all the other 
elements) and from the P-invcntory as determined for [Fc/H] = —3. The parameter /£tl(E) is the fraction of an element E in a star or in the P-inventory 
contributed by SNe II(P) and II(L). The VMS yield pattern is represented by logevMs(E) as calculated from /jji(E) and the observed loge(E) in a star or 
logep(E). For convenience of comparison, the logeyMS values have been shifted to pass through the same Fe abundance of logevMs(F e ) = loget(Fe) = 5.03. 

a Atomic number. 
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